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Abstract 

In "Classical Electrodynamics" |T] a theorem is proved on the average 
of an electrostatic or magnetostatic field over a spherical volume. The 
proof of the theorem is based on an expansion in spherical harmonics 
and it is useful for deriving the singular behaviour of dipole fields. In 
this paper we give a simple proof of this theorem, the key element being 
to use reciprocity. The method also highlights the relation between the 
dipole singularity and the depolarizing field of a sphere, which is thus also 
found very easily. In addition this link enables us to extend the theorem 
to a more general ellipsoidal volume. Unfortunately one can only prove 
the existence of the depolarizing tensor in this way and one still needs 
another method for actually calculating the tensor. Finally the theory is 
also extended to an anisotropic background medium. 

1 Introduction. 

In [T] it is proved that the average of the electrostatic field taken over a spher- 
ical volume equals the electric field at the center of the sphere if the charge 
distribution causing the field is completely located outside the sphere. If on the 
other hand the charge distribution is completely enclosed by the sphere then the 
average is proportional to the elecric dipole moment of the charge distribution 
calculated with respect to the center of the sphere. The proof of this property is 



(a) The charge is located outside of the 
sphere. 



(b) The charge is located inside 
the sphere. 



Figure 1: The average over a spherical volume of the electrostatic field of 
a monopole q is equal to the electrostatic field of a uniform charge density 
p = —q/\Cl\ in the sphere and measured at the location of the monopole. 

based on an expansion in spherical harmonics of the potential of a monopole and 
similar results hold for magnetostatic fields. The second form of the theorem 
has then be used to derive the singularity of the electric/magnetic dipole fields, 
which has received some attention recently [2j[3]. 

In |2]|3]we will show that both theorems can be proved rather easily by invoking 
the reciprocity of the fields and especially for the electrostatic case (^2j the 
derivation is almost immediate. In ^3]the corresponding magnetostatic theorem 
is proved using the conventional current loop model for a magnetic dipole. In 
Q we briefly recall how the singularities of the dipole fields can be found using 
these theorems. In |5]the close link is revealed between these singularities and 
the well-known depolarizing field of a sphere. This also enables us to extend the 
results to more general shapes of the volume and finally in ^6] to an anisotropic 
background medium. 

2 Electrostatic field. 

It is of course sufficient to prove the theorem for the field of a single monopole 
q. Using superposition the extension to an arbitrary charge distribution is then 
straightforward. We thus consider the average of the field of a monopole q, 
located at P, over a spherical volume f2 (see figure [TJ. Using reciprocity we 
conclude that the average field equals the field E p (P) at the source 

point P generated by a charge — q uniformly distributed over the volume with 
density p = —q/\Cl\, where |f2| is the size of this volume. This first step in the 
reasoning is valid for an arbitrary volume. For a spherical volume with center at 
Q the latter field is easily calculated. If the source point P is outside the sphere 
(P ^ il) then E p (P) equals the field of a single monopole — q at the center Q 
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of the sphere and applying reciprocity again it follows that the averaged field 
equals the field of the original charge q measured at the center of the sphere: 

(E q ) n = E p (P) = E_ q (P) = E q (Q) P^n (1) 

The two outer equalities are due to reciprocity, the central one is due to the 
spherical symmetry. 

If on the other hand the source point is enclosed by the sphere (P 6 n) then 
invoking again the spherical symmetry a standard application of Gauss integral 
law yields: 

where p q is the dipole moment of q with respect to the center of the sphere. 



3 Magnetostatic field. 

A similar reasoning can be followed for the magnetostatic case, although the 
calculations are slightly more complicated if the conventional electric-current 
model is used for the magnetic dipole (m = ia with i the current in the loop 
and a the oriented area of the loop). The magnetostatic result can be obtained 
much more economically by using the magnetic-charge model |4] instead (see {j3|. 
However since the current model is more widely adopted and for completeness 
we show the derivation for this model too. 

We consider a <5-like current element which we denote by K (with the current 
density J(f) = KS(f — rp)). Without loosing generality we choose the origin 
at the center Q of the sphere and the z-axis along the direction of K with 
K = Kl z (see figure [2}. Using the reciprocity of Biot-Savart's law the magnetic 
induction (B r )q of the current element at P averaged over the sphere ft equals 
the induction Bj(P) measured at P of a uniform current density J — — K/\n\ 
present in the sphere. The solution of the latter problem is analogous with the 
electrostatic problem with the correspondence q/eo <-> [IqK and E r B v / sin# 
where E r = —dtfi/dr is the radial electrostatic field and B v = — sm8dA z /dr the 
tangential component of the induction, </> being the electrostatic potential and 
A = A z l z the vector potential, and where is the angle between the z-axis and 
the position vector r. It then follows readily that: 

(B R ) n = Bj(P)=B_ it (P) = B R (Q) Pin (3) 

for the situation in figure |5Ja) and for the situation in figure [5Jb) we find simi- 
larly: 

<**>„ ^ B J{ P) = -^IPQI sin^ = -i^PQ x K = ^ P e fl 

(4) 
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(a) The current element is located outside the 
sphere. 



(b) The current element is lo- 
cated inside the sphere. 



Figure 2: The average over a spherical volume fl of the magnetostatic field of 
a current element K is equal to the magnetostatic field of a uniform current 
density J = —K/\Cl\ in the sphere and measured at the location of the current 
element. 



where is the magnetic moment of the current element with respect to the 
center Q of the sphere. Note that the extra factor 2 compared with the electrical 
dipole is due to the factor 1/2 in the definition of the magnetic moment, this 
being a second order moment, whereas the electric dipole moment is of first 
order. The results (|3|_and ([1} are again easily extended for an arbitrary current 
density distribution J(r). 



4 Singularities of the dipole fields. 

Applying the second form of the theorem Q to the field of a dipole p enclosed by 
an arbitrary sphere it follows that the field must contain a singularity — ^r8(f). 
To take this singularity into account the field is then written explicitly as [TJ: 



- p 3rr-r 2 I p 

£ ° E p=^- — 70 ( 5 ) 



where I is the unit tensor and where we have written the field in a form which 
is more amenable to a generalization to anisotropic media (see Sjfj}. The inter- 
pretation of this expression requires some care. If the sphere is centered on the 
dipole then the average of the first non-singular term vanishes only if one ex- 
cludes an infinitesimal spherical volume around the dipole. This actually means 
that (|5j is only valid if the dipole is modelled as an infinitesimally small po- 
larized sphere, as will become more explicit in |5] As an additional bonus we 
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find that, since the average of the field is, according to independent of the 
position of the dipole within the sphere, the principal value integral of the first 
term must also vanish if the sphere is not centered on the dipole. 

The corresponding magnetostatic result can be obtained immediately by using 
the magnetic-charge model H] for the magnetic dipole. In this model a magnetic 
moment is composed of 2 magnetic charges ±q* a distance d apart. Then 
/Ltom = q*d and due to the analogy with the electrostatic case the magnetic field 
H will show a singularity — y<5(f). Since B/fiQ = H + M and M — mS(f) it 
follows that B/fio must have a singularity |m5(r) which agrees with |4]) found 
higher. Finally we note that for the electrostatic case the electric induction D 
also shows a singularity |p5(r). 



5 Depolarizing fields. 

Applying the theorem in (JlJ once more to a dipol^jp we find that the field of 
a dipole averaged over a sphere (not enclosing the dipole) is equal to the field 
of the same dipole measured at the center of the sphere. Applying reciprocity 
to both sides of this equality it follows that the (external) field of a uniformly 
polarized sphere (P — p/|fi|) is equal to the field of a dipole p at its center. The 
same reasoning can be applied for any higher order multipole. Similarly applying 
([2]) for a dipole but now enclosed by the sphere it follows that the electric field 
within a uniformly polarized sphere (P = p/\Q\) is uniform and equal to — 
This is the well-known depolarizing field of a sphere which is thus shown to be 
directly related to the singularity in §5§ and this leads us to the conclusion that 
([5]) is only valid if the dipole is modelled as an infinitesimally small polarized 
sphere. Applying the same reasoning to any higher-order multipole and since 
the dipole moment of any higher-order multipole vanishes we conclude that in 
a sphere filled with a uniform higher-order multipole density the internal field 
vanishes. 

These results can be formalized as follows. Using the same notation as before 
the field of a multipole averaged over an arbitrary volume can be found by 
filling with a corresponding multipole density and looking for the field at the 
source point P: 

(E q ) n = E p (P) 

(E p ) n = Ep(P) (6) 



1 In what follows we limit ourselves to the electrostatic case, but all results are readily 
extended to the magnetic case. 
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where p = — g/|f2|, P = p/|fi|, Q = — g/|f2|, ... are the charge density, polariza- 
tion density and quadrupole density, ... corresponding with the charge q, the 
dipole p, the quadrupole q ... . The multipole fields on the left sides follow from 
each other by differentiation: 

E p =-p- VE q /q 

E= = n-VVE q /q (7) 



Due to reciprocity similar relations hold for the fields on the right sides: 

Ep = -P- X7E p /p 

E= = 5 • VV£ p /p (8) 



with V operating on the position vector of the (original) source point P. Note 
that Q and |8| are valid for arbitrary volumes. The field of a uniformly polar- 
ized body in particular can thus be derived from the field of the same uniformly 
charged body. This method for calculating the field of a uniformly polarized (or 
magnetized) volume has been attributed to Poisson j5]- 

For a sphere in addition the symmetry can be invoked and the external field 
of a uniformly charged sphere simply equals that of a monopole. Therefore the 
field of a uniformly polarized sphere also equals a dipole field and similar results 
hold for higher order multipole moments. Since the internal field of a uniformly 
charged sphere is linear in the position, as shown in ([2}, the internal field of 
a uniformly polarized sphere is then according to the first equation in (JsJ , a 
constant and given by: 

^ = -3^ = -^=^ P ^ (9) 



where we stress once again the relation between the singularity of the dipole 
field on the left and the depolarizing field of a uniformly polarized sphere on the 
right. Due to the linear dependence on the position vector it follows immediately 
that the internal field vanishes for all higher-order uniform multipole densities. 

Having found this relation it becomes straightforward to look for extensions 
of the theorems ([TJ and (|2j to other than spherical volumes. It is obvious 
that an extension to ellipsoidal volumes must be possible for an internal point. 
Therefore we consider firstly a (centric) hollow sphere. From the linear position 
dependence in ([2J it is clear that the averaged electric field of any charge enclosed 
by the hollow sphere is zero and this is equivalent with the well-known property 
that the field of a uniformly charged hollow sphere is zero in the enclosed space. 
It also readily follows that this remains true for an enclosed dipole distribution 



6 



(and for the corresponding polarized hollow sphere) even if the inner sphere is 
eccentric. This argument can be reversed as follows. If the electric field in a 
uniformly charged (centric) hollow sphere vanishes within the enclosed space 
then the internal electric field of a uniformly charged (single) sphere must be 
linear dependent on the position. Indeed if the latter field (say of the outer 
sphere) is denoted by E(r), with the origin taken at the center, then in the 
enclosed space of a hollow sphere the field equals E(r) — sE(r/s), where s is the 
scale factor scaling the outer sphere to the inner one. If this field vanishes then 
E(f) must be linear in r. 

Next we transform a spherical shell into a more general ellipsoidal shell by 
directional scaling. Taking into account the quadratic dependence of the electric 
field on distance it is easy to show that also in the cavity of a uniformly charged 
hollow ellipsoicj^] the field vanishes [6] . Therefore applying the same scaling 
argument as above, the internal field of a uniformly charged ellipsoid is also 
linear in the position (referred to the center) and ([2} can then be written for an 
arbitrary ellipsoid as: 



(E g ) u = E P (P) = ■ QP = -L i{ A) • ^ Pen (10) 

Where Li depends on the shape of the ellipsoid, defined by A : r¥ = 1, A being 
a symmetric positive definite tensor, and the index i refers to "isotropic" since 
we are still considering an isotropic background. Using ^ it follows that: 

(E p ) n = Ep(P) = -%(A) Pen (11) 

showing that Li is the usual depolarizing tensor. Of course one still needs a 
method for calculating this tensor. Today the best known approach is probably 
by the method "separation of variables" in ellipsoidal coordinates j7], but Li 
can also be calculated directly [5]. Perhaps the most elegant expression can be 
found by solving the problem in Fourier-space [9] HQ] : 

%(A) = ~ ■ [ J^dn (12) 

4lT J A:rm 



where n is a unit vector in the (spatial) frequency space and the integration is 
over the surface of the unit sphere. The expression immediately reveals that the 
trace of the depolarizing tensor equals 1. 



Now turning again to a hollow ellipsoidal shell it follows from (11 1 that the field 



is still zero if the inner ellipsoid becomes eccentric, but still aligned with the 



2 By a hollow ellipsoid we mean that inner and outer bounding surfaces are similar ellipsoids, 
which should not be confused with the case that the ellipsoids are confocal. 
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similar outer ellipsoid. This means that if such an eccentric hollow ellipsoid 
encloses a dipole, the average of the field over its volume vanishes. In particular 
the principal value integral of the dipole field over an ellipsoidal volume thus 
vanishes if the exclusion volume is a uniformly scaled version of the ellipsoid 
and irrespective of the position of the dipol^] The corresponding singularity of 
the dipole field is then — L ■ ^S(r) and more in general (JsJ is written as: 



47r r° 

In this case the dipole is thus modelled more generally as an infinitesimally 
small ellipsoid. Obviously there is no unique singularity for a dipole field and 
the singularity depends on the symmetry attributed to the dipole. 

The "external" theorem (jlj cannot be extended from a spherical body to an el- 
lipsoidal body since the field outside a uniformly charged ellipsoid is not longer 
equal to that of a monopole, because the even order multipole moments do not 
vanish. Consequently the field outside a uniformly polarized ellipsoid is also not 
longer equal to a pure dipole field. However if the ellipsoid is made infinitesi- 
mally small, preserving the shape and keeping the dipole moment constant, then 
all higher order moments of the uniform charge distribution will vanish and the 
external field of the polarized body once more becomes a pure dipole field, inde- 
pendent of the shape of the ellipsoid. This explains why the non-singular part 



in (13 1 does not depend on the shape. 



6 Anisotropic medium. 

The results can also be extended to a uniform anisotropic background medium. 
E.g. consider a uniform dielectric medium with symmetrical relative dielectric 
tensor e. The extension is easily handled by introducing a suitable coordinate 
transformation which transforms the anisotropic medium into an isotropic one 
and then applying the results for the isotropic medium [12]. If we preserve the 
electric potential and charge density in corresponding points then the relevant 

3 In an appendix of [TT] this is proved directly and from this the constancy of Li is derived. 
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transformation formulas are given by: 



—1/2 

r = T ■ r 

— =-1/2 

E = T ■ E 

— =1/2 _/ 

D = T D 

- =1/2 =1/2 

e = T ■ f ■ T 

p = |f|l/2f V2 .^ 

— — -1/2 —I — -1/2 

A = T -A ■ T 

— —1/2 —I =-1/2 

L = T -L - T (14) 



where |T| is the determinant of T. With T = e the transformed dielectric tensor 
becomes isotropic (e = /) . Consider now as volume Q an ellipsoid isomorphous 
with the "index ellipsoid" of the medium e : rf = 1, then is transformed into 
a sphere fi' in an isotropic medium, to which ([TJ can be applied. Transforming 
back we find that ([lj remains valid for an anisotropic medium when applied to a 

volume isomorphous with the "index ellipsoid". The depolarizing tensor L(e, A) 
depends on the anisotropy of the medium and on the shape of the ellipsoidal 
volume and follows directly from the last two transformation formula's in ([141: 



l(f,I) = f 1/2 .l l (f 1/2 .5.f 1/2 ).r 1/2 (15) 



and (10 1 and (111 can thus be extended to: 



e Q e-(E q ) n = -L(e,A)-^r Pe!J (16) 



P 
\Q\ 

showing that the second theorem actually holds for the electric displacement 



Eo i-(%> Q = -i(i,3)-4 



instead of for the electric field. From ( 15 I it also follows that the trace of 
L(f, A) is still 1. 

As a special case we consider now a sphere (A — I) and we then obtain in 



_ —1/2 = — 1/2 — 

particular L(e, I) —e ■ Li(e) -e — Li(e) as given by (12 I and where in the 
last step we used the property that Li(e) and e share the same principal axes. It 
follows that e and A separately have the same effect on the depolarizing tensor. 
We consider then an alternative transformation which transforms the ellipsoid 
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into a sphere (T = A and A — I) and obtain in this way an alternative 
expression for the depolarizing tensor: 



= _ = =-1/2 = =1/2 _ =1/2 =1/2 

L(e,A) = A -Li(A -e-A ) • A 



(17) 



The expression (17) with (12) has been further extended to a bianisotropic 
medium 1131. 



Comparing (15) and (17 1 it follows that: 



L(e,A) ■ e = e ■ L(A,e) 

3 •!(?,!) = t(A,t) -3 



(18) 



showing the more general symmetry between anisotropy (e) and shape (A). 
The non-singular part of the dipole field can also be obtained with the first 
coordinate transformation resulting into the following dipole field: 



47r|e| 1 /2 



(=-x _\ 5 / 2 f=-i _\ 
I s : rr I I e : rr I 



3/2 



L(e,A)-p8{¥) (19) 



which reduces to ( 13 1 for e = I. Usually this is called a Green's dyadic and the 



field of an arbitrary polarization density P(r) is then obtained by a convolution 
with this dyadic, with the last term resulting into —L(e, A) ■ P which is only 
present if the observation point is in the source region. The link between this 
contribution and the depolarizing tensor has been noted by Weiglhofer |14] . 
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